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The Landau theory of phase transitions has been re-examined under the framework of a modified 
mean field theory in ferroelectrics. By doing so, one can see that there are two atomic movements 
involved in the ferroelectric phase transition; the first corresponds to the vibration of crystalline 
lattice, which will render phonon mode softening at the critical point, and the second represents 
the slow evolution of a partially ordered nematic phase formed by the cooperative behavior of high- 
temperature structure precursors. In this hierarchical dynamic structure, the former fast dynamics 
could be significantly modulated by the latter slow dynamics in the vicinity of the Curie temperature; 
it then turns out that it is the behavior of the nematic phase on approaching the critical point that 
makes the Landau theory deviate from experimental observations. 
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In 1937, Landau formulated an elegant theory, in which 
he proposed the concept of broken symmetry to explain 
continuous phase transitions; Landau also introduced a 
unique thermodynamic variable, i.e., the order param- 
eter, to demonstrate his broken symmetry concept 
For the phase transition from paramagnetism to ferro- 
magnctism, the spontaneous magnetization M is defined 
as the order parameter. Thus, within the framework of 
the Landau theory, M has the following values: M = if 
temperature T is above the critical temperature T c and 
M = ± constant if T is below T c ; the former (T > T c ) 
corresponds to one phase with higher symmetry and the 
latter (T < T c ) another phase with lower symmetry. In 
the absence of the external magnetic field and in the 
vicinity of T c , the Landau free energy F per unit vol- 
ume can be written as follows [l[ 
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where Fq is the free energy that is independent of M; 
both a and b are coefficients; a is also defined as a = 
a t = a T ^ Tc , in which t = T ~ Tc is called the reduced 
temperature and ag is a coefficient. The Landau theory 
is one of the basic theoretical tools for describing phase 
transitions. By using the broken symmetry concept and 
the Landau free energy, continuous phase transitions can 
be explained qualitatively and quantitatively and the uni- 
versal behavior of certain critical exponents near T c can 
be defined. For the ferromagnetic phase transition, the 
following critical exponents are given by the Landau the- 
ory: 
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here C and Xm represent specific heat and magnetic sus- 
ceptibility, respectively; H m represents the external mag- 
netic field. Unfortunately, these estimated results are 
not in good agreement with experimental observations. 
The corresponding experimental data and the exponents 
given by the Landau theory are summarized in Table 1 
for convenience. 

The major criticism of the Landau theory can be con- 
cisely summarized as follows: thermal fluctuations are 
neglected in the Landau theory; however, near the criti- 
cal temperature, thermal fluctuations become important 
and render the Landau theory inaccurate 0. Before we 
discuss this common belief, it is perhaps worth briefly 
discussing the mean field approximation of the Landau 
theory so that one can see what else, in addition to ther- 
mal fluctuations, is missing in the Landau theory. 

Let us consider a uniaxial magnetic material with the 
simple cubic structure. Its Hamiltonian can be written, 
via the Ising model, as 
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where H and Sj = ±1 represent the Hamiltonian and 
the Ising spin on the zth lattice site, respectively; Sj is 
defined as the neighboring Ising spin of the ith site; h is 
the external magnetic field on the ith site; J represents 
the interaction strength between the nearest neighbors. 
If we assume that the mean value of spin, m, is site- 
independent, i.e., to = (sj) = (sj), then 5si = Si — m 
and Ssj = Sj — m represent the random values caused 
by thermal fluctuations. Substitute both Sj and Sj by 
Sj = to + Ssi and Sj = m + Ssj in the above equation, 
the Hamiltonian can be simplified as follows Q 
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(2) and the corresponding partition function Z can also be 



BJ /J 

where N is the total number of lattice sites; z is the 
number of nearest neighbors of each spin; ks represents 
the Boltzmann constant. Using Eqs. @ and (JSJ) , one 
can further determine that m = when T > T c or 
m = ± constant when T <T C Q. After such a simplifi- 
cation, one can see that the Hamiltonian of a many-body 
system given in Eq. can be reduced to that of an 
effective one-body system given in Eq. in which the 
second order fluctuation term, SsiSsj, is neglected and 
Ssi and 5s j are split and integrated into h e ff = h + zJm; 
here, h e ff represents an effective mean field or molecular 
field on the ith site, which arises both from the exter- 
nal field, h, and from the exchange field induced by the 
neighboring spins, zJm. The expression of h e ff can be 
further extended to the corresponding effective macro- 
scopic mean field, H e tf = H m + aM; here a is a coef- 
ficient and M represents the total of m within the ma- 
terial. The above mean field approximation of the Ising 
model reflects a part of the nature of the Landau theory, 
which is based on an assumption that the considered ma- 
terial undergoes harmonic vibrations when perturbed by 
external stimuli (including thermal fluctuations). Under 
this situation, the atoms sitting on lattice sites or the 
corresponding Ising spins on lattice sites of the consid- 
ered material can be regarded as distinguishable physical 
quantities that obey the Boltzmann distribution; the mu- 
tual interactions between the atoms or the spins due to 
thermal fluctuations can be regarded as random variables 
and their net effect can be neglected. Thus, the above- 
mentioned simplified Hamiltonian and the partition func- 
tion can be derived. Within this mean field approxima- 
tion, the Landau theory assumes that the net value of 
total spins is a site-independent non-zero constant value 
M and the mutual interactions between spins are zero 
when T < T c or both of them are zero when T > T c . 

However, the lattice vibration of the considered ma- 
terial is by no means harmonic all the time, especially 
when approaching the critical point. If the anharmonic 
lattice vibration induced by thermal fluctuations makes 
the crystalline lattice distorted locally in the considered 
material, the assumption of the Landau theory is still 
valid. In this case, the spins on the distorted lattice sites 
might not obey the Boltzmann distribution but the net 
value of those spins and the mutual interactions between 
them can be assumed to be zero since those distorted 
lattice sites are randomly distributed in the spatial do- 
main. Within the framework of the Landau theory, the 
corresponding Hamiltonian has a global symmetry since 
it is supposed to be invariant with respect to spatially 
uniform group operations. For those spins on the dis- 
torted lattice sites and the mutual interactions between 
them, their net contribution to the magnetization is zero 
and thus they do not affect the global symmetry of the 
Hamiltonian. If, for some reason, the spins and the in- 
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FIG. 1. Diagrammatic representation of the emergence of 
distorted crystal lattice; (a) the undeformed magnetic crystal 
lattice; (b) and (c) the distorted magnetic crystal lattice; (d) 
the undeformed crystal lattice of a perovskite material in the 
paraelectric state, in which there is no net dipole moment; (e) 
and (f ) the distorted perovskite crystal lattice with net dipole 
moment represented by arrows. 



teractions cannot be averaged out and do form a global 
magnetization in both T < T c and T > T c cases (of 
course, such a magnetization, if exists, must be weak), 
they can indeed affect the global symmetry and then fur- 
ther shift the phase transition behavior of the considered 
material. This situation has not been considered in the 
Landau theory. Furthermore, there is another concern in 
the Landau theory. In thermodynamics, material prop- 
erties are described by using continuous functions so that 
the considered material must be treated as a continuum 
Q. Since the Landau theory is a thermodynamic model, 
it must be subject to this thermodynamic limit; such 
an essential prerequisite is implicitly assumed to be only 
missing at T c but satisfied above and below T c in the 
Landau theory due to the assumption of the global sym- 
metry. It later becomes clear that it is these two factors 
that render the Landau theory inaccurate and alter the 
asymptotic behavior of the corresponding physical quan- 
tities defined in Eqs. fl2J|. 

Low temperature Nematic phase Isotropic 

phase phase 



FIG. 2. Schematic representation of the cooperative behav- 
ior spectrum of HTSPs; below temperature Tn, the quantity 
of disordered structures is too small to behave cooperatively; 
above temperature Tni, thermal energy is large enough to de- 
stroy any possible cooperative behavior of disordered struc- 
tures; between Tn and Tni, a nematic phase as the ensemble 
of disordered structures exists. 

Since the above-mentioned spins and interactions are 
related to the distorted crystal lattice, we need to take a 
close look at its formation and the possible dynamic be- 
havior of the corresponding microscopic disordered struc- 
tures under thermal fluctuations. For a crystalline or 
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polycrystallinc material, it has a purely ordered state at 
temperatures just above absolute zero and a completely 
disordered state above its melting point; at temperatures 
between those two extremes, the material should have 
both ordered and disordered states and the formation of 
the latter is largely due to thermal fluctuations (for sim- 
plicity, we here do not consider defects or microscopic 
disordered structures generated during material manu- 
facturing processes). For the normal crystal lattice of 
magnetic materials shown in Fig. 1 (a), even at tem- 
peratures far below the melting point, there always ex- 
ists the probability that certain atoms could gain extra 
kinetic energy from thermal fluctuations to move quasi- 
permanently away from their original equilibrium posi- 
tions and make the crystal lattice distorted as shown in 
Fig. 1 (b) or (c). Such distorted lattice could change 
the local electron degeneracy and then makes the local 
spins and the effective local magnetic field very compli- 
cated in the considered magnetic material. In order to 
simplify the theoretical analysis to be discussed, we will 
mainly consider perovskitc materials and the ferroelec- 
tric phase transition in the rest of the paper. For sim- 
plicity, we will first consider a perovskitc material in the 
paraclectric state. For the undeformed crystal lattice of 
such a material shown in Fig. 1 (d), one can imagine 
that the centers of gravity of the negative and the posi- 
tive ions coincide leading to zero net dipole moment; for 
the distorted lattice shown in Fig. 1 (e) or (f), how- 
ever, there is relative displacement between the centers 
due to the local non-uniform deformation, which results 
in net dipole moment. The local disordered structures 
corresponding to such distorted crystal lattice are de- 
fined as high-temperature structure precursors (HTSPs) . 
Since the formation of the distorted structures is caused 
by thermal fluctuations, their quantity is proportional to 
temperature. Therefore, as temperature rises, the quan- 
tity also increases. At a certain temperature TV, the 
quantity of HTSPs has reached a threshold or the effec- 
tive distance between HTSPs has been reduced below a 
critical value so that HTSPs will start to interact with 
each other, which can be regarded as their cooperative 
or self-organization behavior, and then form a unique 
nematic phase. The driving force behind this structural 
transformation is the competition between energy and 
entropy. Let us consider a crystalline perovskite mate- 
rial containing HTSPs; its Gibbs free energy in the ab- 
sence of external fields can be written as G = U — TS. 
Clearly, the internal energy U will increase due to the in- 
crement of the local strain and electric potential energy 
generated by HTSPs shown in Fig. 1 (e) or (f). If T re- 
mains unchanged or changes slowly, the entropy S must 
rise to reduce G. The simplest way to increase S is that 
the chemical bonds between the atoms on the distorted 
lattice sites are partially broken so that the correspond- 
ing disordered structures could gain more freedom to ro- 
tate and are oriented along local preferred directions to 
form a nematic phase, which leads to a decrease in the 
orientational entropy but an increase in the positional 



entropy and, eventually, results in a net increase in the 
total entropy. In addition, since the nematic phase has 
no polarity (the cylindrical symmetry Dooh) [6[, the in- 
duced electric potential energy shown in Fig. 1 (e) or (f) 
will decrease when the phase is formed so that G could 
be further reduced. If we use v to represent the local 
preferred direction and approximately treat HTSPs as 
molecules, then their orientational order parameter, S op , 
can be written as [H, @ 

S op = i(3 (v\n) 2 I) = i((3cos 2 ^ - l)>, (6) 

where ( } represents the average; v l is defined as the 
given direction of the disordered structure located at the 
position i; ft is usually called the director that represents 
a particular direction; 6 l is defined as the angle between 
v l and n at the position i. If < S op < I, we can say 
that HTSPs cooperatively form a nematic phase. S op = 1 
corresponds to an ideal case, in which all HTSPs are 
perfectly aligned. If temperature continues to rise, at a 
certain point T = Tjyj, the thermal energy will be large 
enough to disturb HTSPs, which makes HTSPs randomly 
oriented and then renders S op = 0. Therefore, when T > 
T/v/, all HTSPs behave like a normal liquid, which has 
an isotropic phase. The cooperative behavior spectrum 
of HTSPs is shown schematically in Fig. 2. 
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FIG. 3. Schematic representation of how the dynamic be- 
havior of k could alter phase transitions; (a) k — > 0.5~ co- 
incides with Tc in an allowed temperature range, then the 
ferroelectricity-to-paraelectricity phase transition will occur 
at Tc; (b) k — > 0.5~ occurs far before T reaches Tc, then 
there exists no second order phase transitions. 



Wc now consider the behavior of this nematic phase in 
a perovskitc material in the ferroelectric state. Since HT- 
SPs are usually metastable, the built-in polarization field 
Eb (in the absence of the external electric field) would 
tend to disturb them. According to Le Chatelier's prin- 
ciple, the nematic phase would then undergo a specific 
structural change to counteract any imposed change by 
the filed Q. Therefore, the effective built-in polarization 
field inside the specimen will be E e J^ = Ej, — gPy L = 
gP — gPh', here the negative sign is due to Le Chatelier's 
principle; both g and g are coefficients, respectively; both 
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P and Ph represent the spontaneous polarization (the or- 
der parameter) and the induced counteractive polariza- 
tion of the nematic phase, respectively. For the sake of 
convenience, we can further write E^f* = g(l — k)P and 
define k as 

k Q S op (T — T N ) 
k = ^ {T N <T < Tat/), (7) 

-iJV 

where both k and kg are dimcnsionlcss coefficients. For a 
unit-volume ferroelectric material containing HTSPs, in 
the absence of the external electric field, the correspond- 
ing free energy F up to the fourth order of P can be 
written as 

F = Fa + \a{l - k)P 2 - ^ Vh P^ + h(l - k) 2 P 4 , (8) 

rjh is a coefficient. This equation represents a special case 
of the Landau-de Gennes free energy Q , which indirectly 
demonstrates the co-existence of the ferroelectric phase 
transition and a weak first order structural transforma- 
tion. If r)h is very small and can be neglected, then the 
above equation will reduce to a modified Landau free en- 
ergy, which is written below 

F = F a + ia(l - k)P 2 + h(l - fc) 2 P 4 

„ „ / ak o b, , bk „,\ 
= F +F L + (- Y P + ~ A k P - Y P ) ' (9) 

where Fl = f-P 2 + \P A ■ In the above equation, one can 
see that k actually represents, in the statistical sense, the 
fraction of the electric potential energy, which is gener- 
ated by the cooperative movement of HTSPs to reduce 
the potential energy brought by the order parameter. 
Now we can briefly summarize what we have done so 
far: the potential energy corresponding to both the crys- 
talline phase and the nematic phase as the ensemble of 
HTSPs or, more precisely, the mutual interactions be- 
tween the atoms of HTSPs is integrated together by using 
a modified mean field theory (MMFT), in which HTSPs 
are treated as a partially ordered liquid. In this way, two 
phases, the crystalline phase and the nematic phase, must 
be considered together since the latter is constrained by 
the former. The dynamics of the crystalline phase is fast, 
which corresponds to the vibration of crystalline lattice. 
On the other hand, the symmetry group in the nematic 
phase is the continuous rotation group, which has an un- 
countable continuum of symmetry elements that renders 
the corresponding relaxation extremely long and, thus, 
its dynamic behavior is extremely slow Q. Therefore, 
the integration of these two phases does not violate Lan- 
dau's global symmetry assumption in ferroelectric mate- 
rials and his concept of broken symmetry in the ferro- 
electric phase transition; however, it does embed a weak 
first order structural transformation in the second order 
phase transition, which is partially reflected by the term 
(1 — k) in the modified Landau free energy given in Eq. 
© • Before discussing how this term could alter the phase 



transition behavior, we need to find out the meaning of 
k. For a crystalline ferroelectric material, k = when 
T < Tjy, which corresponds to the situation that Eq. 
© will reduce to the regular Landau free energy, and 
k will increase when T rises. When Tjy < T < T^j, 
the nematic phase of HTSPs will start to emerge in the 
considered material and the chemical bonds between cer- 
tain atoms of HTSPs are partially broken so that they 
cannot be considered the integral part of the crystalline 
lattice of the considered material. Therefore, k = 0.5 
represents, in the statistical sense, the thermodynamic 
limit; below this limit (k < 0.5), the material is still a 
continuum; above this limit (k > 0.5), the material can- 
not be regarded as a continuum anymore Thus, from 
the viewpoint of thermodynamics, Eq. Q is only valid 
when k < 0.5. Such a situation has been reflected in the 
Coffin-Manson law, which is written below 

| « e f (2N f )~^ = e f (2N f )-^, (10) 

where e p is the induced plastic strain; e/ is defined as 
the ductility coefficient; Nf represents the fatigue life, 
and /3cm, which has been proved to be $cm = 1 — k 
is called the Coffin-Manson exponent. It has been 
observed that Pcm possesses the remarkable universal- 
ity; Pcm ~ 0.5 has been found in single-phased metallic 
materials (see Ref. and the references cited therein). 
Pcm > 0.5 or k < 0.5 can be interpreted that the test 
specimen undergoes elastic deformation (or the specimen 
is still a continuum); whereas (3cm < 0.5 or k > 0.5 
demonstrates that fracture and crack have emerged and 
started to grow in the specimen (in other words, the spec- 
imen is not a continuum anymore). Therefore, the Coffin- 
Manson law actually describes the fatigue behavior near 
the limit where fracture and crack will emerge in the 
specimen. 

We now consider the situation, k — > 0.5, in the 
ferroelectricity-to-paraelcctricity phase transition. From 
the above discussion, it is clear that the crystalline phase 
and the nematic phase form a hierarchical structure. At 
very low temperatures, the atomic movement of the for- 
mer is frozen, then the latter does not exist; as tem- 
perature rises or the movement of the former becomes 
fast, the latter starts to emerge. Therefore, the latter is 
constrained by the former. At the critical point T c , the 
crystalline phase collapses due to phonon mode softening 
and the considered ferroelectric material must be in the 
disordered state. The nematic phase could quickly grow 
up if the temperature does not rise too fast. Thus, k 
must reach the thermodynamic limit at some point near 
T c (without loss of generality, we assume < T c < T^i 
as shown in Fig. 2). This situation brings a rather subtle 
change to the phase transition. As shown in Fig. 3, if 
k — > 0.5~ coincides with Tc in a very small temperature 
range, then the ferroelectricity-to-paraelectricity phase 
transition will occur; if, however, k — > 0.5~ occurs far 
before T reaches T^, then there exists no second order 
phase transitions. We shall demonstrate how the dy- 
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namic behavior of k would alter the phase transition and 
the corresponding critical exponents in the former case. 

By taking the partial derivative of F given in Eq. ([9]) 
with respect to P and then letting the result be equal to 
zero, we have the following equation 



dF 
dP 



(1 - k)P [a + 6(1 - k)P 2 } = 0. 



(11) 



One of the solutions of this equation is given below: 

t 



P = ± 



b 



1 - k 



t < (12) 



Contrary to the Landau theory, the phase transition char- 
acterized by the free energy given in Eq. (|9]) involves two 
atomic movements near T c as previously mentioned; the 
result of the fast dynamic behavior can be interpreted 
by the soft-mode theory [13, [LH but it could be signif- 
icantly modulated by slowly fluctuating HTSPs in the 
nematic phase in the vicinity of T c . Therefore, near 
T c , the critical exponent (3 is determined by both move- 
ments as shown in the above equation, in which both t 
and 1 — k have key influence on /?. To estimate j3, one 
has to do the following modification. We first rewrite 



i-fc 
l+fc 



P 1 
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(1 - k)P 2 as (1 - k)P 2 = 1 + (1 + k) 

near T c , P <C 1, k — > 1 (mathematically, k approaches 
1 at T c but, physically, k — > 0.5 when T is reaching 
Tc; beyond this limit, the equations written here can- 
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not be used), and jj^ < 1, thus 

taking advantage of the binomial series expansion, we 
have the following relationship: (1 — k)P 2 w P^ +2fe ; here 

P l = T^ P \-Th + l = U p2 + Tk- E ^ ©can 
then be rewritten as 



p m = ±(-^) 2+2t M)^. t<o 

The mathematical expression of j3 is given by 

1 



2 + 2/fc- 



(13) 
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where k = k — > 0.5 . If k = 0, (3 will become fihandau- 
Thus, the upper limit of (3 is ^; since k < 0.5, its lower 
limit should be i. 

At the critical point (t = 0), we apply a weak electric 
field E to the considered material, the exact differential 
of its Gibbs free energy will be dF = -SdT + EdD = 
—SdT+EdP; here D represents the electric displacement 
and D = e$E + P, eo is the electric permittivity of free 
space. E is then written as 



E = \JP) T = 5(1 _ k)2p3 ™ bP ™ ik - (15) 
Thus, the critical exponent 8 can be defined as 



S = 3 + 4k. 



(16) 



If we let k = 0, then S will reduce to 5L anc [ au . The lower 
limit of 6 could then be regarded as 3; since k < 0.5, its 
upper limit should be 5. 

By taking the partial derivative of E with respect to 
P, we get the following formula 



dPj T ~ \dP 2 ), 



a(l-k)+3b(l-k) 2 P A . (17) 



Then the corresponding electric susceptibility Xe is given 
below 



1 dP 1 



1 



Xe = — 



s dE e a{l-k) + 3b(l-k) 2 P 2 ' 



(18) 



When t > 0, the order parameter P will be zero. How- 
ever, the polarization generated by the nematic phase will 
not be zero. In the vicinity of T c , k — > 5; thus the nematic 
phase has grown up, due to Le Chatelier's principle, and 
its volume is too large to respond to the applied weak 
electric field so that the remaining counteractive polar- 
ization docs not contribute to \e and can be neglected. 
Then Eq. (fTg]) will reduce to 
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When t < 0, by using P 2 
modified as 



h(l-fc) ; 



(19) 



Eq. (|T8l) can be 
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Since a <C 1 in the vicinity of T c , we can modify a(l — k) 
by exploiting the method previously used. We let a(l — 



*) = ! + (! + k) 
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By using the binomial 



series expansion, we have the following relationship o(l ■ 



k) w aj+ k ; here a m = |-p|< 
flT9l is then rewritten as 
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The critical exponent 7 can be defined as 
7 = l + k+, 



. Eq. 



(21) 



(22) 



where k + = k — > 0.5 + . Similarly, we can also derive 
the expression of the critical exponent 7', which is given 
below 



7' = l + k~ 



(23) 



If k = 0, both 7 and 7' will reduce to jLandau and 
^Landaui respectively; thus, the lower limit of both 7 and 
7' is 1. If we let k = 0.5, we can get their upper limit, 
which is equal to 1.5. In addition, k + > fc~, thus 7 is 
slightly larger than 7'. 

We now consider the asymptotic behavior of the spe- 
cific heat C, which is defined as C ~ — T§j§-, in the 
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TABLE I. Critical exponents of continuous phase transitions; the values with the superscript symbol [a] in the following table 
are taken from Table 3.1 of Ref. 0| and the references cited therein; the ones with [b] are taken from Table 5.4.2 of Ref. [3] 
and the references cited therein; the ones with [c] are the simulation results of the three-dimensional renormalization-group of 
Ising systems reported in Refs. 0, [f| . Here MMFT and RG represent the modified mean field theory proposed in this research 
paper and the renormalization group, respectively. 

Exponents Experiment RG (n = 1) Landau theory MMFT (k = 0.47) Upper & lower limits 

a 0.11~0.12 N 0.110 [cl 0.113 l>a>0 

P 0.28 ~ 0.38 [i,) 0.325 [c] ± 0.34 ± > /3 > ± 
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2 2 ^ r " 3 

1.20 ~ 1.47 [tl 1.241 [c] 1 1.47 1.5 > 7 > 1 



8 4.6 ~ 4.9 W 4.82 w 3 4.88 5 > S > 3 



vicinity of T c . We first rewrite C as: ^- = — if-f 1 



-t 



d 2 F 



g T 2 



gag d 2 F 
' T? da 2 



If k = 0, by using F 







g T 2 Q T 2 ~ T 2 Q a 2 T 2 Q a 2 

given in Eq. ([9]), we can calculate C as follows: |f 

whcn t > and jr = ^26T°" wnen * < 0. Clearly, 

in both cases, C does not diverge and only has a discon- 

a 2 

tinuity jbW a t ^ = 0. Thus, the corresponding critical 
exponents are usually defined as a Landau = a' Landau = 0. 

If k ^ 0, however, both a and a' will not be zero due 
to the embedded first order structural transformation. 
When t > 0, the order parameter P given in Eq. © will 
be zero; however, the counteractive polarization gener- 
ated by the nematic phase will not be zero. By replacing 
P 2 with — gfpzfcl , we can rewrite the free energy based 
on the remaining counteractive polarization as follows 



F = F 



where h 



2(l-fc) ' 



l+k„ 



1 



b 4(1 - k) 2 
Let ah = a(l 



F 



1 



(ah) 2 

k m ) = 1 
2(l-fc) 



(24) 



k 



(1 + 

; near T c , 



a«l and k — > 1 as explained previously, thus 2 ^ fc ^ 



and 



2(l-fc) 
k 



< 1. Similarly, using the binomial 



series expansion, we can get the following relationship: 



ah 



a 



l+k„ 



3fc-2 



and a m = a — 



l , here fc m 1 2(1— k) 

h 1 = a + k~ 2 ' Then the above equation can be 

rewritten as 



F^F 



a 



2+2k„ 



Ct2 Tp Cj2 T 

It can be seen that = 5-4 
can be calculated as 

aa Q d 2 F a§ d 2 F 



(25) 



thus the specific heat C 



C 



aa d 2 F 

T c da 2 ~ Y c ~diJ = 'l^daj, 

a m ' 



a 2 d 2 F 
T c dal 



where J\ 



a (2 + 2k^(l+2k m ) ^ 



3k-2 
k 



(26) 



Now 



we have two parameters describing the slow and the fast 



asymptotic behavior of C, respectively; the former is 1 + 
2k m and the latter 2k m . In practice, we usually measure 
the slow asymptotic behavior; therefore, we can define 
the critical exponent a as 



a = -(l + 2k m ) 



l-2fc+ 
1 - k+ ' 



(27) 



When t < 0, the corresponding free energy can be written 
as 



F = F + F L 



(-a) 2 



k' 2 



b 4(1 -k) 2 



(28) 



In this case, Fl does not affect the critical behavior (see 
a Landau ~ 0) an< ^ can ^ e neglected. Thus, the critical 
exponent a' could be derived, in the same manner, as 



a 



1 - 2k- 
l-k~ 



(29) 



If k = 0, the values of both a and a' will become a 
constant value of 1, which is their upper limit; if we let 
k = 0.5, we will get their lower limit, which is equal to 
0. Wc now try to determine which exponent, a or a', 
is larger. Subtracting a' from a, we have the following 
expression 



k+ 



a — a 



(i-k-)(i-k+y 



(30) 



Since fc~ < fc + , therefore, a' should be slightly larger 
than a. 

So far, the mathematical expressions of the critical ex- 
ponents have been derived from the ferroelectric phase 
transition. But, by using the same method, they can 
also be derived from other continuous phase transitions. 
We now try to estimate the exact values of the expo- 
nents. In the vicinity of T c , the value of k of a perfect 
single crystal should be slightly less than 0.5, which is its 
theoretical limit. In practice, due to the existence of de- 
fects, k varies from material to material near the critical 
point. For most magnetic materials with relative simple 
crystalline structures, k should be close to 0.5. However, 
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for some perovskite systems with complicated crystalline 
structures, k may vary over a wide range of values. For 
instance, a Spain-Ukraine group recently investigated the 
evolution of the ferroelectric transition with Se doping in 
a ferroelectric family Sn2P2(Se :E Si_ 2 ;)6 (0 < x < 1); for 
x varying from 0.2 to 0.3, the exponent a they measured 
ranges from 0.14 to 0.34 [l2j, which corresponds to the 
values of k varying from 0.4624 to 0.3976. In this pa- 
per, without loss of generality, we let k « 0.47; then the 
corresponding critical exponents can be calculated and 
the results are summarized in Table 1. Clearly, our re- 
sults are in good agreement with the experimental data, 
which were mainly obtained in magnetic materials, and 
the results of renormalization group calculations. 

It might be interesting to briefly discuss the Ginzburg 
criterion [l3|, [H[ > m which the spatial correlation of the 
order parameter is used to quantitatively determine when 
the Landau theory is valid. The problem with this crite- 
rion is that, near T c , the dynamic behavior of the corre- 
lation could be significantly influenced and modulated by 
the cooperative movement of HTSPs, which is not con- 
sidered in Ginzburg's model. Therefore, the Ginzburg 
criterion may not provide a complete physical picture on 
why the Landau theory is not valid under certain circum- 
stances. 

Concluding remarks - the theoretical consideration and 
derivation presented in this paper have shown that there 
are two atomic movements involved in the ferroelectric 
phase transition. These two dynamics correspond to the 
vibration of crystalline lattice and the evolution of a par- 
tially ordered nematic phase formed by the cooperative 
behavior of high-temperature structure precursors, re- 
spectively; they are not independent of each other but 
rather constitute a hierarchical dynamic structure. The 



vibration represents the fast dynamics and the evolution 
corresponds to the slow dynamics, which is constrained 
by the vibration. The slow dynamics not only "imposes" 
the thermodynamic limit on the occurrence of broken 
symmetry near the critical point but also alters the phase 
transition behavior, which are demonstrated by the fact 
that the cooperative behavior of HTSPs could certainly 
shift the asymptotic behavior of the corresponding physi- 
cal quantities defined in Eqs. J2]). Now, it is perhaps safe 
to say that it is the behavior of the slow dynamics that 
makes the Landau theory deviate from experimental ob- 
servations. The studies presented in this paper also rise 
a question whether it is necessary to draw a strict sep- 
arating line between first-order and second-order phase 
transitions. According to the classical definition, there 
is no latent heat involved in second order phase tran- 
sitions 15|; within the framework of the theory devel- 
oped in this paper, however, the slow dynamics always 
involves a weak first order structural transformation so 
that, contrary to the Landau theory, there should be la- 
tent heat involved in continuous phase transitions. Fi- 
nally, we would like to emphasize that the conclusions 
given in this paper should apply to general continuous 
phase transitions though they are drawn from the ferro- 
electric phase transition. 
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